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Abstract 

The Yudin-Ivanov formula [Phys. Rev. A 64, 013409 (2001)] is generalized such that the most 
general analytical expression for single-electron spectra, which includes the dependence on the 
instantaneous laser phase, is obtained within the strong field approximation. No assumptions on 
the momentum of the electron are made. Previously known formulas for single-electron spectra 
can be obtained as approximations to the general formula. 
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Although an initial theoretical understanding of strong field ionization was put forth by 
Keldysh [1] as early as 1964, many questions remain unresolved. As far as single-electron 
ionization in the presence of a linearly polarized laser field is concerned, there are two 
important topics. The first, namely the ionization rate as a function of instantaneous laser 
phase, was studied in depth by Yudin and Ivanov [2] (see also [3, 4]). However, their 
result assumes zero initial momentum of the liberated electron. Effects due to nonzero 
initial momentum have yet to be included. The second topic pertains to the single-electron 
spectra, that is, the ionization rate as a function of the final momentum of an electron. 
Despite much work on this topic (see, for example, Refs. [5-7] and references therein) a 
universally accepted formula is lacking and the discussion is still ongoing. Among the most 
accurate results, Goreslavskii et al. [6] have obtained an expression for the complete single- 
electron ionization spectrum, but without consideration of laser phase. The present paper 
derives a more general formula [see Eq. (8)] that includes the dependences on both the 
instantaneous laser phase and the final electron momentum. 

It is very convenient to formulate the Keldysh theory in terms of the Dykhne adiabatic 
approximation, as was presented for in [8]. According to the Dykhne method [9-11] (see 
also [8, 12]), if the Hamiltonian of a system H(t) is a slowly varying function of time t, and 
H{t)ipn{'t) — En{t)ipn{'t) , {n — i, f), then the probability F of the transition ipi ipf is given 
by (atomic units h — nie — \e\ — 1 are used throughout) 



where ti is any point on the real axis of t, and to is the transition point, i.e., a complex root 
of the equation 



which lies in the upper half-plane. If there are several roots, we must choose one that is the 
closest to the real axis of t. Moreover, there are no assumptions regarding the form of the 
Hamiltonian. 

Now we shall apply the Dykhne approach to the problem of ionization of a single electron 
under the infiuence of a linearly polarized laser field with the frequency cu and the strength 
F. The initial and final energies for such a process are given by 




(1) 



Ei{to) — Ef{to), 



(2) 



Ei{t) 



(3a) 
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E,(t) = l[k + A(t)f, 



(3b) 



where Ip is the ionization potential, k is the canonical momentum (measured on the detec- 
tor), and A{t) — — ^sin(a;t). 

According to Eq. (1), the probability of one-electron ionization F can be written as 



r = exp [-2lvciS{k\\MJp)] , 
where S is the action. Equation (2) can be rewritten in terms of 5" as 



(4) 



(5) 



Note that the analogy between the saddle point S-matrix calculations [13], where transitions 
are calculated using stationary points of the action, and the Dykhne approach can be seen 
from Eqs. (4) and (5). The transition point is given by 



utn = Arcsin 



where 7 is the Keldysh parameter 
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7 = jV^I, - 



2C/„ 



and Up — (^)^ is the ponderomotive potential. In order to extract the imaginary and real 
parts of this solution, the following equation [14] can be used 



Arcsin {x + iy) — 2Kt^ -\- arcsin (3 



a 



+ Vq;2 - 1 



(7) 



where K is an integer and 

{^} = ^\^(^+W7± ^\/(^^W7- 

Using Eq. (7) in Eq. (4), we obtain 

r(7, /cy, /cj_) cx exp 



-^/(7,fc||,M 



(8) 
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where 



f{^,kii,k±) = ^1 + 2^ + arccosh a - Va^ - 1 + 



g [1 - 2/3^ 
2^ 




a 

A;2 = k\\^ + k^_\ 

Note that « > 1. It must be stressed here that no assumptions on the momentum of the 
electron have been made. However, Eq. (8) has an exponential accuracy because the influence 
of the Coulomb field of a nucleus cannot be accounted for by the strong field approximation. 
The correct exponential prefactor has been obtained within the Perelomov-Popov-Terent'ev 
(PPT) approach [15-18]. 

Similarly to the Yudin-Ivanov formula, Eq. (8) is vahd if the strength of the laser field F 
depends on time, F — > Eog{t), where the envelope g{t) of the pulse is assumed to be nearly 
constant during one-half of a laser cycle. 

Equation (8) is the central result of this paper. In the following, this equation is applied 
to some special cases in order to establish connections with previous results. 

In the case of zero final momentum (A; = 0), we have that a — a/I + 7^ and (3 — 0. In 
this hmit we recover the original Keldysh formula [1] 



/(7, 0, 0) = ( 1 + ^i^") arcsinhT - ^^^^^ 
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In the tunneling limit (7^1) the following formulas can be obtained. Expanding the 
function /(7, k\\, k±) in a Taylor series up to third order with respect to 7 and setting k± = 0, 
we obtain 



^11 3 

3a; ^ 



(9) 



r(7,A;||,0)«r(7,0,0)exp 

Eq. (9) has been derived by a classical approach in Ref. [19] (see also Ref. [20]). Discussions 
regarding the physical origin of Eq. (9) are presented in Ref. [13]. Performing the same 
expansion and setting k\\ = 0, we obtain 

2 (^^^ + 2/,)=^/^ 



r(7, 0, k±) oc exp 



(10) 
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This equation has been derived in Ref. [20]. For small values of k±, Eq. (10) can be 
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approximated by 



r(7,0,A;x) Rir(7,0,0)exp 



(11) 



Let us fix k± — and continue working in the tunneling regime. For the case of high 
kinetic energy, such as k\\'^/2 > 2Up and ^ /{AUp) — 1 » 7, we obtain 



/9~1, 



and the ionization rate F is given by 



F (X exp ^ 



2C/„ 



+ 1 I arccosh 



4C/„ 



-3i 



- 1 



(12) 



Eq. (12) has been obtained in Ref. [5]. 

Calculating the asymptotic expansion of the function 7(7, A;||,0) for k\\'/2 ^ 2Up, we 
obtain 

'F2exp(3)\ 2- 



Foe 



(A;||72 > 2Up, k± = 0). 



(13) 



Eq. (13) has been obtained for tunneling ionization in Ref. [5]. Here, we have proved that 
Eq. (13) is valid for arbitrary 7. A similar formula can be derived for k± ^ 



Foe 



{kj_y2 > 2Up, k\\ = 0), 



(14) 



which is also valid for arbitrary values of the Keldysh parameter 7. 

Consider the asymptotic expansion of Eq. (8) for large values of Ip {Ip ^ k'^/2). In this 
case a and (3 can be approximated by 

7 k\\ 



a -x/l + 7^, 
Using these equations, we obtain 



k' 



^1 + 7^^ 



/(7, ^||, ^±) ~ /(7, 0, 0) + — arcsinh7 
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k\\ 



(15) 



\/l + 7^ 2/p 

Eq. (15) has been reached within the PPT approach [15-18] (see also [7]). 

As mentioned above, Goreslavskii et al. [6] have derived an expression for the spectral- 
angular distribution of single-electron ionization without any assumptions on the momentum 
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of the electron. However, they have summed over saddlepoints, i.e., the contribution from 
previous laser cycles has been taken into account. On the contrary, we have not performed 
any summation because we are interested in the most recent contribution to ionization. 

Therefore, our result does account for the phase dependence of the ionization rate, unlike 
thatofRef. [6]. 

To make the phase dependance explicit in Eq. (8), we apply the substitution 

kii^kii-A{t). (16) 

The analytical expression for the ionization rate as a function of a laser phase when A: = 
has been achieved by Yuding and Ivanov [2]. Thus, Eq. (8) is seen to be a generalization of 
the Yudin-Ivanov formula. 

Note that generally speaking, there is no unique and consistent way of defining the in- 
stantaneous ionization rates within quantum mechanics, and such a definition is a topic of an 
ongoing discussion (see, e.g., Refs. [21, 22] and references therein). However, the instanta- 
neous ionization rates are indeed rigorously defined within the quasiclassical approximation 
(the Yudin-Ivanov formula), and we have employed this approach in the current paper. Al- 
ternatively, one can approximate the instantaneous ionization rates by the static ionization 
rates at each point in time using the instantaneous value of the laser field. 

Lastly, we illustrate Eq. (8) for the case of a hydrogen atom. The single-electron ioniza- 
tion spectra in the multiphoton regime (7 ^ 1) and in the tunneling regime (7 <^ 1) are 
plotted in Fig. l.A and Fig. l.B respectively. One concludes that the smaller 7, the more 
elongated the single-electron spectrum. We can notice that the maxima of both the spectra 
are at the origin. Nevertheless, a dip at the origin has been observed experimentally [23, 24] 
in the parallel-momentum distribution for the nobel gases within the tunneling regime, and 
afterwards it has been investigated theoretically in Ref. [25] and references therein. How- 
ever, such a phenomenon is beyond Eq. (8). The phase dependence of ionization for different 
initial momenta, recovered by means of Eq. (16), is illustrated in Fig. 2 for selected positive 
momenta. The curves for negative momenta are mirror reflections (through the axis = 0) 
of the corresponding positive curves. Figure 3 shows that the cutoff of the single-electron 
spectrum in the tunneling regime (the dashed line) corresponds exactly to the kinetic en- 
ergy 2Up, which is the maximum kinetic energy of a classical electron oscillating under the 
influence of a linearly polarized laser field. 
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FIG. 1: The single-electron spectrum of a hydrogen atom at 800 nm; plot A: at 1 x 10^^ W/cm^ 
(7 = 3.376); plot B at 6 x 10^^ W/cm^ (7 = 0.4357). 

In summary, we have derived a novel expression for strong field ionization including both 
the dependence on momenta and instantaneous laser phase. Previous results [1, 2, 5, 20] 
regarding strong field ionization can be recovered as special cases of the general formula. The 
present result concerns only the exponential dependence of the ionization process. However, 
supplementing our approach with the pre-exponential factor taken from PPT [7, 15-18], the 
present result is the most general formula obtained within the quasi-lassical approximation. 

The author is grateful to M. Spanner, G. L. Yudin, and M. Yu. Ivanov for many valuable 
discussions and remarks, and for outstanding comments on the manuscript. 
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FIG. 2: The plot of r(7, k + ^ sin</), 0) for a hydrogen atom at 800 nm; plot A: at 1 x lO^'^ W/cm' 
(7 = 3.376); plot B at 6 x 10^^ W/cm^ (7 = 0.4357). 
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